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Abstract 

We compute the attenuation of gamma-rays from re- 
mote quasars due to interactions with ambient low-energy 
photons produced by stars in galaxies and show that the 
next-generation gamma ray telescopes can probe the cos- 
mic star formation rate at high redshifts. 

Key words: galaxies: evolution - galaxies: diffuse background 
radiation - quasars: TeV gamma-rays 



1. Introduction 

The cosmic star formation rate (CSFR) at high redshifts 
is a key probe to determine the time-table of structure 
formation in the universe and hence to further constrain 
cosmology. Number counts of galaxies only give access to 
small parts of the sky and are subject to severe selec- 
tion effects depending on the wavelength range of the ob- 
servational campaign such as, for example, extinction by 
dust using deep photometry in the optical. Furthermore, 
absolute sky photometry is biased due to galactic fore- 
grounds making it difficult to assess truly diffuse extra- 
galactic background light (EBL) which could arise due to 
the decay of unstable relics from the early universe. There- 
fore, an independent method unaffected by these obstacles 
is required, and we suggest along with Salamon & Stecker 
(1998) and Primack et al. (1999) to use the attenuation of 
gamma rays from distant quasars for this purpose. 

EGRET has detected that many flat-spectrum radio 
quasars at high redshifts emit gamma rays at least up 
to 10 GeV showing generally no indication for a spectral 
break. A few nearby flat-spectrum radio sources have even 
been detected up to TeV (10 12 eV) energies using ground- 
based Cherenkov telescopes suggesting that at least some 
of the sources exhibit smooth power-law spectra in the 
GeV-to-TeV energy range. The gamma rays of energy E 1 
can interact with ambient low-energy photons of energy e 
creating electron-positron pairs provided the threshold en- 
ergy E y ~ (m e c 2 ) 2 /e is exceeded. As a result, the gamma 
ray sources should show a cutoff in their spectrum at the 
energy where the mean free path for pair creation be- 
comes shorter than the distance to the sources, i.e. where 
r 77 (i? 7 ) ~ cr 77 ri c bi[(TO e c 2 ) 2 /£' 7 ]d z > 1 where a 77 ~ ctt/3, 
n e bi denotes the photon number density of the EBL, and 



d z is the cosmological distance of the gamma ray source. 
Generally, since the number density of background pho- 
tons increases with decreasing photon energy, and thus 
increasing gamma ray energy above threshold, the more 
distant gamma ray sources must show lower cutoff ener- 
gies than the closer gamma ray sources and the precise 
running of the cutoff energies with redshift will tell us 
about the shape and evolution of the EBL. 

Our strategy now is to develop a simple model for 
the redshift-dependent EBL based on the observed (and 
extinction-corrected) CSFR and then to determine the 
corresponding values of gamma ray cutoff energies. We 
start with a description of the galaxy evolution templates 
used to determine the EBL when folded with the CSFR, 
and proceed with the computation of the gamma ray op- 
tical depth. Finally, we discuss that possible deviations 
of future measurements of the gamma ray cutoff ener- 
gies would require changes in the EBL model (e.g. in the 
CSFR) or force us to consider additional sources to con- 
tribute to the EBL. 



2. Spectral energy distribution 
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Figure 1. Spectral energy distribution of a simple stellar 
population with mass 1O 11 M0 and solar metalicity evolv- 
ing from an age of 0.01 Gyrs to 17 Gyrs. 
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We use the spectral synthesis model calculated by Moller 
et al. (1999), which describes the evolution of the spec- 
tral energy distribution (SED) of a single burst popula- 
tion (SSP) with time (see Figure Q). The stars have solar 
metalicity and Salpeter IMF. We simplify the temporal 
evolution by reducing it to only two spectral states. The 
first state lasts from ignition to an age of 0.28 Gyrs. We 
call this state the early phase and it is dominated by the 
ultraviolet and optical radiation of young massive stars. 
The second state is defined to begin at an age of 0.28 Gyrs. 
We call this state the late phase and it is dominated by 
the optical and infrared light of the older stars in the pop- 
ulation. To get the spectrum representative for each phase 
we calculate the corresponding mean SED from the tem- 
plates. Note that the late phase lasts over a cosmologically 
significant portion of time and therefore must be weighted 
with the CSFR. Furthermore, we include the effect of the 
interstellar medium (ISM) by assuming an isotropic and 
homogeneous gas and dust distribution around the stars of 
the SSP. We use Osterbrook's Case B for optical thick gas 
clouds, i.e. total absorbtion of photons below the wave- 
length of 911 A and reemission of the absorbed power in 
the optical regime by bremsstrahlung and H a line emis- 
sion. The dust is treated the same way as in Desert et 
al. (1990) considering a three component model with big 
grains, small grains and PAH's. For the absorbtion we use 
an extinction curve similar to the one of the SMC. The 
reemission is described by four Planck spectra with two 
temperatures for each phase reprensenting warm and cold 
dust components. The resulting spectra including absorb- 
tion and re-emission are shown in figure 0. 
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Figure 2. Spectral energy distributions of a simple stellar 
population for the early (the first 0.28 Gyrs) and late (af- 
ter 0.28 Gyrs) phases. The solid line includes (dashed line 
excludes ) absorbtion and reemission of the ISM. 



3. Cosmic star formation rate 

Now let us look at galaxy evolution on global scales. In our 
model, the stars in galaxies are all represented by SSPs. So 
far we did not include a particular star formation rate, and 
we will do this now by using the cosmic star formation rate 
(CSFR) as determined from galaxy counts in the optical 
(e.g., Madau 1996). Since dust and gas absorbs optical 
photons, the original Madau CSFR underestimates the 
true SFR in particular at high redshifts by at least a factor 
of ~ 3 (Pettini et al. 1998) or maybe even ~ 4.7 (Steidel et 
al. 1999). All data agree, however, in a CSFR increasing 
up to a redshift z — ~ 1 — 2 and decreasing at z > Zb, 
see Fig. ||. We adopt in the following a power-law fit to 
this CSFR shown as the solid line in Fig. || 

p\{z) oc (l + z) a (1) 

with a = a m > for z < and a = f3 m < for z > z&. 

1.000 i i 




0.0 1.0 2.0 3.0 4.0 5.0 6.0 7.0 



FigureS. Cosmic star formation rate P* vs. redshift z. The 
data are taken from Madau (1996; solid circles), corrected 
by Pettini et al. (1998; open squares), Gallego (1995; dia- 
mond) and further data from Steidel et al. (1999; stars). 
The solid line indicates the fit used in our model 



4. Emissivity 
Solving the convolution integral 

£ x (t)= f F x (t-t')p* (t')dt' 
Jo 

we can calculate the emissivity at different cosmic times 
t = t(z). For the comparison with measurments we will 
integrate over redshift z and change to the comoving emis- 
sivity 

£ x (z) = (1 + z)- 2 / T x (t{z) - t{z')) p t {z')—dz'. (2) 
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We used the parameters in Tab. 1 to obtain the results 
shown in Fig. The emissivity £\(z) is plotted for three 
different wavelengths demonstrating satisfactory agreement 
Note that in particular the steeper increase at 0.28 /im and 
the more shallow increase at 1.0 /im are well reproduced 
by the model. 

E%_ v E l B _y dm Pm Z b ^[Mq^' 1 Mpc" 3 ]" 

0.1 0.03 3.4 -2 1.3 0.13 
Table 1. Choosen parameter for the calculation of the 
emissivity. Note that we used the better value [3 m = — 1 
for the CSFR (Fig. 3) and the background (Fig. 5) calcu- 
lation. 
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Figure 4-- Emissivity vs. redshift . Data taken from Lilly et 
al. (1996; squares), Conolly et al. (1997; circles) and Ellis 
et al. (1996; dots). 



5. EXTRAGALACTIC BACKGROUND LIGHT 

The intensity of the EBL follows from the emissivity 

xh = ,±jy y(z) !?L(±±l}\, (3) 

We compute the background intensity at several red- 
shifts and show the results in Fig. IE. Beyond a redshift 
of z ~ 1.3, the EBL decreases because of the decreas- 
ing CSFR and the expansion of the universe. The shape 
of the spectrum at high redshifts is determined by the 
early phase spectrum, whereas for lower redshifts the late 
phase becomes more and more important. This is a conse- 
quence of of the maximum star formation rate at redshifts 



z — 1 — 2 from which the stellar populations have consid- 
erably aged until the present. There are no direct mea- 
surements of the EBL at high redshift (with the exception 
of inferences of the UV background based on the proxim- 
ity effect which are consistent with our model predictions 
if we employ empirical power law extensions of the UV 
spectra beyond the Lyman edge). We can, however, com- 
pare our results with the present-day data obtaining good 
agreement as shown in Fig. |6|. 
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Figure 5. Diffuse Extragalactic background light due to 
stars in galaxies, for various redshifts. 



6. Gamma ray attenuation 

The gamma-ray attenuation can be parametrized by the 
optical depth for photon-photon pair production 

T{R 1 ,z q )=cj^ J o jf —t nchx {e,z'){l + zf- (4) 

■a 11 {E 1 , e, /Lt, z')de d\i dz 

Using Eq. (3) for the intensity of the background pho- 
tons we can calculate the number density n e bi(e, z'). In- 
serting n e bi(e, z') in Eq. (4) we can compute the opti- 
cal depth for various redshifts. The optical depth enters 
the attenuation factor exp[— t(Ej, z q )] for gamma rays 
(see Fig. [?]) modifying quasar power law spectra. Future 
gamma ray observations with ground-based air Cherenkov 
telescopes in the energy range 10 GeV to 10 TeV can ef- 
ficiently probe the run of the CSFR with redshift. E.g., 
the faster the decrease above the maximum of the CSFR 
(the larger the slope (3 m in Eq. (1)), the less attenuation 
occurs for quasars at redshifts z»l On the other hand, 
a plateauish CSFR combined with the effect of dust keeps 
the UV radiation density low at high redshifts, but leads 
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to a strong FIR bump in the EBL which leads to strong 
attenuation of > 10 TeV gamma rays from nearby z<l 
radio sources. 
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Figure!. Attenuation factor exp[— t{E 1i z q )\ for gamma- 
rays as a function of gamma ray energy E 7 for various 
quasar redshifts z q . 



Figure 6. Comparison of the model with the observed 
present-day EBL. Data points refer to Armand et al. 
(1994; open triangle), Pozetti et al. (1998; squares /HST), 
Dwek et al. (1998; filled circle at 35/im/DIRBE), Altieri 
et al. (1998; triangle /ISO CAM), and Hauser et al. (1998; 
diamond/IRAS, blue filled circles/ DIRBE). 



7. Summary and conclusions 

In order to facilitate a comparison of gamma ray cutoff 
energies in the spectra of high-redshift quasars with the- 
oretical predictions of gamma ray attenuation due to pair 
creation, we developed a simple model for the evolving 
EBL based on the observed CSFR taking into account ex- 
tinction and re-emission by dust and gas. Further model 
inputs are the IMF, the average metalicity, relative contri- 
butions of warm and cold dust to the emission of early- and 
late-type galaxies, and the Hubble constant for which we 
have adopted values in close agreement with observations. 
The CSFR at high redshifts is not accurately determined 
by the data, and we have tried to use a value which is con- 
sistent with both the present-day EBL and high-redshift 
galaxy counts. So far we have neither considered the ex- 
tension of the galaxy spectra beyond the Lyman edge, the 
metalicity dependence of galaxy spectra, nor the contri- 
bution of quasars to the EBL. The predicted attenuation 
of gamma rays from nearby radio sources such as Mrk501 
or Mrk421 leads to a quasi-exponential cutoff near 10 TeV 
due to interactions with the infrared part of the EBL, and 
to a cutoff near 100 GeV due to interactions with the op- 
tical/UV part of the EBL for sources with redshifts in 
excess of one. 
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